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We present a general interpolation theory for the phenomenological effects of thermal fluctuations 
in superconductors. Fluctuations are described by a simple gauge invariant extension of the gaussian 
effective potential for the Ginzburg-Landau static model. The approach is shown to be a genuine 
variational method, and to be stationary for infinitesimal gauge variations around the Landau gauge. 
Correlation and penetration lengths are shown to depart from the mean field behaviour in a more 
or less wide range of temperature below the critical regime, depending on the class of material 
considered. The method is quite general and yields a very good interpolation of the experimental 
data for very different materials. 

PACS numbers: 74.20.De, 74.25. Bt, 11.15.Tk 

I. INTRODUCTION 

In a recent paperQ we have shown that the Gaussian Effective Potential (GEP) can describe the behaviour of 
superconductors, thus aUowing for a comparison with the experimental data. The GEP has a long historyf 



superco 



^ JTJ and has been discussed by several authors as a tool for describing the breaking of symmetry in the 

standard model of electroweak interactions 0, ^3 1 but in that context no real comparison with experimental data 
will be achievable until the detection of the Higgs boson. Thus High Tc superconductors represent the best way to 
test the reliability of the method. While the comparison with experimental data is of special importance as a test for 
the GEP variational method itself, the method can be regarded as a general interpolation scheme for the pre-critical 
region of superconductors: too close to the critical point for thermal fluctuations being negligible, but not yet in the 
critical range where universality sets in. In this pre-critical range, the GEP provides an interpolation between the 
mean- field behaviour and the critical limit. Besides standard phenomenological parameters the method relies on one 
free parameter which characterizes the width of the pre-critical range and allows for a very good fit of the experimental 
data for very different materials ranging from cuprates like YBa2Cu207_5 to the double band superconductor MgB2. 
The free parameter turns out to be a characteristic energy scale of the material. 

The general phenomenology of superconductivity can be described by the standard Anderson-Higgs mechanism: the 
supercurrent is carried by pairs of charged fermions whose non- vanishing expectation value breaks the gauge symmetry, 
and endows the gauge bosons with a non-zero mass. Thus, the standard Ginzburg-Landau (GL) effective Lagrangian 
provides the best framework for a general description of the phenomenology. Moreover, as the GL action can be seen 
as a power expansion of the exact action around the critical point, the GL action must be recovered by any microscopic 
theory at least around the transition. Therefore, whatever the origin of the microscopic pairing mechanism, the GL 
action is a sound starting point for a general interpolation scheme of superconductivity. Of course we cannot trust the 
mean-field approach to the GL effective theory, and we expect that in high Tc superconductors many unconventional 
properties should have to do with the breaking down of the simple mean-field picture. Actually, unconventional 
superconductors are characterized by a very small correlation length ^, which allows the experimentalists to get closer 
^ to the critical point where the thermal fiuctuations cannot be neglected and the mean-field approximation is doomed 
k>( know there is no clear evidence that the critical universal behaviour can always be reached in 

rN the cuprate superconductors: a weakly charged superfluid regime has been predicted and observed^ ITsl IT^ in 
' nearly optimally doped samples, while a charged critical behaviour has been reported for the underdoped samples |l^. 
■ ■ ■ ' Recently the superconducting pyrochlores have been shown to open up a window onto the charged critical regime [ig. 

However it is out of doubt that an intermediate range of temperature is now accessible, where thermal fluctuations 
are not negligible even when the sample is still out of the truly critical regime. Thus, in order to describe some 
unconventional properties of high Tc superconductors, we need to incorporate the role of thermal fluctuations, but 
unfortunately we cannot rely on the standard renormalization group methods jl^ which would only describe the 
limiting universal behaviour. We need an approximate approach to the GL action for the non-universal regime 
where the behaviour depends on the physical parameters of the sample, and we would prefer a non-perturbative 
approximation in order to deal with any coupling, whatever its strengthness. 

As a toy model for electroweak interactions, the scalar electrodynamics in four dimensions has been studied by 
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Ibanez-Meier et al.[l2 who computed the GEP by use of Cartesian coordinates for the field components. This choice 
breaks the U(l) symmetry of the field and gives rise to an unphysical, and undesirable, massive degree of freedom. 
However, in three space dimensions, the U(l) scalar electrodynamics is equivalent to the standard static GL effective 
model of superconductivity, and the GEP becomes a variational tool for superconductivity. In this paper, we show 
that the unphysical degree of freedom can be avoided by use of polar coordinates for the field components. Moreover, 
the choice of a U(l) symmetric Gaussian functional restores the gauge symmetry and makes the theory gauge invariant 
for infinitesimal gauge changes around the Landau gauge. In a previous work we discussed the same approximation 
in unitarity gauge(l|. The method has been extended to two-dimensional systems by Abreu et al,20j in order to 
study film superconductors. Unfortunately, we found out that the claimed exact integration of the longitudinal 
component of the electromagnetic vector field is only approximatively exact. Actually, as we show in the present 
paper, that approximation is a variational one and its accuracy can be checked by standard variational arguments. 
Indeed, by a direct comparison of the energies, we show that the present method yields a better result compared with 
the Cartesian coordinate GEP of Ibafiez-Meier et al[T^. 

The variational method provides a way to evaluate important phcnomenological quantities like the correlation 
length ^ and the penetration depth A that emerge as the solution of a set of coupled equations. The GL parameter 
kgl = A/^, whose critical behaviour has been addressed by Herbut and Tesanovic IQj, is here shown to be described 
by a very simple relation in the pre-critical regime where it is found to be temperature-dependent, in contrast to 
the mean- field result. Thus, the GL parameter has been chosen as a useful measure of thermal fluctuations: its 
dependence on temperature is a sign of the breaking down of the mean-field approximation which predicts a constant 
Kgl ■ Unfortunately, experimental data on the temperature dependence of the GL parameter are quite fragmentary 
in the literature. The behaviour of ^ and A can be extracted by the knowledge of the critical magnetic fields, but 
experimental limits on these data have narrowed our chances of a wide comparison with experiments. Up to now, our 
GEP has shown the capability of fitting the GL parameter for all the materials that we have been able to test. 

The paper is organized as follows: in section II the GL action and partition function are reported and the symmetry 
of the model is discussed together with the conditions for the gauge invariance of the variational functional; in section 
III the integration of the longitudinal vector field is shown to be a variational approximation; in section IV the GEP 
is evaluated in polar coordinates and the result is compared with the Cartesian coordinate method of Ref.|l2l|: in 
section V the variational equations are solved for some set of phenomenological parameters and the GL parameter is 
compared with the available experimental data; some comments and final remarks on the interpolation method are 
reported. 

II. THE GL ACTION AND GAUGE INVARIANCE 

Let us consider the standard static GL actionHS 
5 = I d^x 



(1) 



Here is a complex (charged) scalar field, its covariant derivative is defined according to 

D^,=d^+ ieeAf, (2) 

and /i, = 1,2,3 run over the three space dimensions. The components of the magnetic field i^^^j/ = d^Au — d^Af^ 
satisfy 
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^F^,F'^- = |VxA|^ (3) 

and the partition function is defined by the functional integral 

Z = j i?[0,0*,^Je-^. (4) 

The action S in Eq.(Q) has a local U{1) symmetry as it is invariant for a local gauge transformation 

A ^ A + Vx{x) (5) 

^ ^^-leBXix) (6) 
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where xi^) is an arbitrary function. The integration over A is then redundant in Eq.Q and a gauge fixing term 
must be inserted in order to make the partition function finite. According to the standard De Witt-Faddeev-Popov 
method[2^ the partition function can be written as 

Z = I D[^,<j>*,A^]e~''e-'f". (7) 

where the gauge fixing action is 

= /^'^^/' (8) 

and f{A) = is an arbitrary gauge constraint. Z is invariant for any change of the parameter e and of the constraint 
/. With some abuse of language, this invariance property is referred to as gauge invariance while it is a more general 
invariance as Z does not depend on the shape of the weight factor which has been added in Eq.©. Only for e ^ the 
weight factor exp{—Sfix) becomes a J-function which enforces the constraint / = on the vector field A. Thus gauge 
invariance denotes the invariance of the theory for any change of the constraint / = in the limit e — > 0. This is a 
weaker condition, but unfortunately even this is not fulfilled by some approximations to the partition function. In some 
approximation schemes, such as the perturbative method, the action S is split in two parts S = Sq + Si and different 
approximations are considered for the two parts: for instance Sq might need no approximations while exp(— S*!) could 
be expanded as a power series of some small parameter. Even when S is gauge invariant the approximate Z can 
result to be gauge dependent. It is easy to show that if both Sq and Si are gauge invariant then the approximate 
Z is gauge invariant in the sense that in the limit e — ^ any gauge change f f + leaves the approximate Z 
unchanged. Thus in order to obtain a gauge invariant treatment we must take care that the exact action S is split in 
gauge invariant terms. 

The GEP is a sort of optimized first order approximation: in Cartesian coordinates[l2j the action is split according 

to 



Sq ^ / d^x 



(9) 



where ipi + i(/)2 = </> and Si = S — Sq. Here 17,;, A and (pi are arbitrary variational parameters, and the free energy is 
evaluated up to first order in the expansion of Si. The result is then optimized by variation of the free parameters. 
When the shifts ipi are non vanishing the action Sq is not U{1) invariant, and the approximate partition function 
fails to be gauge invariant. The lack of symmetry is evident from the appearance of an unphysical massive Goldstone 
boson (in fact, there are two different mass parmeters fii and 02)- In the next section we show by a variational 
argument that in polar coordinates the action Sq can be taken as 



5*0 = / d^x 



(10) 



where the field is decomposed according to = pexp(i7) and p is a U{1) invariant real field. The explicit integration 
of the phase 7 makes the action invariant for rotations in the (t>i,(t>2 plane. As before fl, A and the shift ip are 
variational parameters. We still have a non-invariant term quadratic in the vector field. However, if we take 



any infinitesimal gauge change f ^ f + Sf yields 
and up to a surface term 



/ = d^Af^ (11) 
A^,^A^ + d^x (12) 



A^A^^A^A,,-2xI + 0{5f) (13) 

so that for e — > the constraint f — restores gauge invariance up to first order in the gauge change around the 
Landau gauge — 0. This discussion motivates in part the choice of polar coordinates that gives rise to an 

explicitly rotational invariant variational functional and no unphysical massive Goldstone boson. Another motivation 
arises from the more pragmatic observation that the free energy is lower with the choice of this polar coordinate 
variational functional, and since we show that both methods are genuine variational methods a lower free energy 
means an higher accuracy of the results. 
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III. GL ACTION IN POLAR COORDINATES 



In the previous section we have given some motivations for the choice of polar coordinates in Landau gauge. In 
three space dimensions that gauge is the transverse gauge / = V ■ A. A simple coordinate change 4> — > pexp(i7) in 
the partition function Eq.Q yields 

d^xC r - / d?xC 



/— I a XL, r — / a xL.^ 
D[A^,p']e J J D[^]e J (14) 

where C is the phase independent lagrangian 

£ = i(V X Af + Id.pd^p+^mlp' + \bP^ + lelp'A^ + 1(V ■ A)^ (15) 
and is the sum of the Lagrangian terms which depend on the phase 7 

A = \p^d^,-id^-i + eBp'^d^.jA''. (16) 

If the last term of £^ were neglected, the phase 7 could be integrated exactly (the integral is Gaussian). It has been 
claimed^, I2(tI l2ll | that in transverse gauge that term should vanish exactly: the argument was that in transverse 
gauge there is no longitudinal component of the vector field. In fact, the Fourier transform of the gradient 9^7 is 
proportional to the wave vector k and this is orthogonal to the Fourier transform of A for the gauge constraint 
/(k) = k • A(k) = 0. However, this holds only for Fourier transforms. The last term of Cy contains a p^ factor which 
is not constant, and thus, in this case, the space integral of this term in the action cannot be replaced by the simple 
inverse-space integral of the product of the Fourier transforms (we could do that if the factor p were constant). Thus 
we cannot neglect that term and the exact integration over the phase 7 would depend on the field A. We stress 
that any gauge change does not solve the problem either, as an undesirable longitudinal component of the field A 
would take the place of the phase 7. For instance in Ref.0| the phase 7 is gauged away but the quadratic term 
e^p^A^ contains an implicit coupling term e^p^A^ • At between the transversal component At and the longitudinal 
component A^ of the vector field. Again the integral of this coupling term would vanish if p were constant (as it is 
obvious by the same Fourier transform argument discussed above), but it is not vanishing in general and prevents 
from an exact integration of the longitudinal component. 

We can get rid of the phase by a variational argument: it turns out that the neglection of the phase terms is not 
an exact integration, but a good variational approximation. Let us denote by £0 the phase dependent Lagrangian 
without the last term 

Co - C^ieB = 0) = ip'a^79^7- (17) 
We observe that up to constant factors the exact integration over 7 yields 

- / d'xCo ___ 1 

\ p 

Thus we may write the Dlp"^] integral in Ea. lfTHl as 



We define the average over the phase as 



/D[7]e-/'''"^' 
With this notation the exact partition function Ea. (|14|) reads 



((•■•))7 = T^TtM- (20) 



d^xC I - / d'^xcBP^dniA'^X 
Z= I D[A^,p\e J (e J ) (21) 
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and the convexity of the exponential function ensures that 



Z > 



D[A^,p]e J e J ^ . 



The average in the right hand side vanishes (it is hnear in 7), and the approximate partition function Zp 



Zp = J D[A^,p]e 



satisfies the variational constraint 

Z > Zp 

so that the approximate free energy JFp = — In Zp is bounded by the exact free energy J- = —InZ 

This bound ensures that, even neglecting the last term of the phase dependent lagrangian in Ea. (|16() . the resulting 
approximate partition function Zp still gives a genuine variational approximation. 



(22) 

(23) 

(24) 
(25) 



IV. THE GEP METHOD IN POLAR COORDINATES 

In this section we study the GEP method for the polar coordinate partition function Eq. (|23|l with the lagrangian C 
defined according to Ea. H15|l . The GEP rnay be evaluated by the same S expansion method discussed in Ref.[l3| and 
[25| and also reported by Camarda et al.|i|. Here the GEP represents a variational estimate of the free energy J^p. 

Inserting a source term for the real field p the partition function Zp reads 



Zp[j] = j D[A^,p]e 
and the free energy is given by the Legendre transform 

J-lf] = — hi Zp + J tfxjip 
where ip is the average value of p. As usual we introduce a shifted field 

and then we split the lagrangian into two parts 



(TxC - j (Pxjp 

e 



(26) 

(27) 
(28) 



(29) 



where Cn is the sum of two free-field terms: a vector field An, with mass A and a real scalar field 6 with mass VL: 



Co 



-(V X Af + -A^A^A^' 



(V-A) 



2e 



The interaction then reads 



+ 1 {eW - A') ApA^ + e\H^ApA^~^ + ^e^A^A^^^ 



where 



(30) 



(31) 

(32) 

(33) 

(34) 

(35) 
(36) 
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The first order expansion of the free energy is 



= iTrln[.g-i(x,2/)] +iTrln[G;J(a;,y)] + 

+ f cPxlvo + V2g{x, x) + 3v4g{x, x)'^ + ie| (g{x, x) + ip'^ - A^) G^,^{x, x) 



(37) 



where g{x,y) is the free-particle propagator for the scalar field, and G^v{x,y) is the free-particle propagator for the 
vector field 



(2^ 



,-ik(x-y) 



5^,,{e + A") + ( i - 1 ) k^k. 



In the limit e — > 0, up to an additive constant 

Trln[G^J(x,2/)] =2V j ^ ln(fc2 + A^) 



(38) 



(39) 



where V is the total volume. Dropping all constant terms, the free energy density 14// = ^/V (effective potential) 
reads 

K//M = /i(r!) + 2/i(A) + 



+ (e|^2^eyo(r!)-A2)/o(A) 
where the divergent integrals /„ are defined according to 

cPk 1 



h{M) 



(27r)3 M2 + A;2 



(40) 



(41) 



(42) 



and are regularized by insertion of a finite cut-off fc < A. 

The free energy H4U|I now depends on the mass parameters $7, A and on the field shift </?. These are the variational 
parameters that must be determined by the minimization of the energy density Vgf / . At the stationary point Vef / 
is the GEP and the mass parameters give the inverse correlation lengths for the fields, the so called coherence length 
^ = l/Sl and penetration depth A — 1/A. The field shift tp is the order parameter of the phase transition: when p ^ 
at the minimum of Ve// the U{1) symmetry is broken in the ground state and the system is in the superconducting 
phase. 

The stationary conditions 



do? 



= 



(43) 



give two coupled gap equations: 



dVe 



eff 



9A2 



(44) 



(45) 



(46) 



For any ip value the equations must be solved numerically, and the minimum point values and A must be inserted 
back into Ea. (|40|) in order to get the Gaussian free energy Veff{p>) as a function of the order parameter p. For a 
negative and small enough m^, we find that Ve/f has a minimum at a non zero value of </? = (pmin > 0, thus indicating 
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that the system is in the broken-symmetry superconducting phase. Of course the masses ft, A only take their physical 
value at the minimum point of the free energy ipmin- This point is found by requiring that 

where as usual the partial derivative is allowed as far as the gap equations H45|) .l|46 |l are satisfied^J- The condition 
(|47|l combined with the gap equation (|45() yields the very simple result 

^Ln = (48) 

However, we notice that here the mass ft must be found by solution of the coupled gap equations. Thus Eas. H48|l . 
(|45|l and 14()|) are a set of coupled equations and must be solved together in order to find the physical values for the 
correlation lengths and the order parameter. 

Insertion of Eq. (|48|l into Eq. H46() yields the simple relation for the GL parameter kq l 

where kq = e^/(8AB) is the mean- field GL parameter which does not depend on temperature. As discussed by 
Camarda et al.|l|, Ea. H49() shows that the GL parameter is predicted to be temperature dependent through the non 
trivial dependence of and ipmin- At low temperature, where the order parameter ifmin is large, the deviation from 
the mean-field value kq is negligible. Conversely, close to the critical point, where the order parameter is vanishing, 
the correction factor in Ea. H49|) becomes very important. The deviation from the mean- field prediction kgl = '^o also 
depends on the cut-off parameter A as for A —> all the effects of fluctuations vanish together with the integrals /„. 
This parameter defines another length scale £ — 1/A, which is a characteristic of the material. 

A test of the present method comes from the direct comparison with the cartesian coordinate GEP discussed by 
Ibanez-Meier et al.|0. Both methods are genuine variational methods and a direct comparison of the free energy 
density gives information on the accuracy. In Fig.l we report our polar coordinate GEP together with the three- 
dimensional version of the cartesian coordinate GEP. The polar coordinate GEP yields a lower free energy density 
thus indicating that, besides being more symmetric and appealing (there are no unphysical degree of freedom and no 
massive gauge bosons), the method is also more effective. 



V. INTERPOLATION OF THE EXPERIMENTAL DATA 



As we mentioned earlier, another appealing aspect of the GEP approximation in three space-time dimensions resides 
in the chance to compare its results with the phenomenology of high Tc superconductors. For this comparison to be 
consistent we need to fix the bare parameters of the GL action; following jj, we can use the standard derivation of 
the GL action (^)(see |2^ ) to find a connection between microscopic first-principle quantities and phenomenologic 
bare parameters. Thus, we have: 



es = T-VkMo (50) 
Tic 

where 2e represents the charge of a Cooper pair, Tc is the critical temperature and the zero temperature coherence 
length. 

The knowledge of the zero temperature coherence length and penetration depth enables us to fix the other param- 
eters; regarding the bare mass parameter to^ as a linear function of temperature, we obtain: 

{ml - ml) (51) 

where TOq is the value which is required in order to find Vl — l/^o from the gap equation (|45|) dX Lp — (pmin, and m^ 
is the value of at the transition point. The mean-field approximation predicts that mj. = 0, but the effect of the 
fluctuations is to shift the transition point to a negative non-vanishing value. 

Finally, the bare coupling is flxed through Ea. H46|) by requiring that at zero temperature (i.e. for = TOq) 
the penetration depth Aq = 1/^- 



m^ 



T 



8 



The cutoff A still remains to be fixed; thus, the interpolation scheme for the superconducting properties we have 
to deal with contains one free parameter (i.e. A), which is a characteristic energy scale of the sample, and will be 
determined by a direct fit of the experimental data. 

We have compared our theoretical results with the experimental GL parameter for three different materials: 
TlsCaaBazCuaOio (T^ = 121.5 K, = 1-36 nm, kq = 100)^2^], YBasCuzOr-a (T^ = 62.2 K, = 1.65nm, 
Ko — 187) jl^l and MgB2 (Tc = 38. 6K, i^o = 3.9nm, kq = 31.66) |28|. For the last two materials we have obtained the 
GL parameter by the relation : 



c2 



2iJci In K 



(52) 



(where Hci and ffc2 represent the lower and upper critical fields, respectively), and extrapolated the zero temperature 
phenomenological values needed to fix the bare parameters by a linear fit. 

The small number of available experimental data is due to the difficulty of performing measurements of the coherence 
length and penetration depth (or, alternatively, of the two critical fields) up to the pre-critical region. 

In Fig. 2, Fig. 3 and Fig. 4 the comparisons between our interpolation curves, as obtained by Ea. H49l) . and the data 
relative to the three materials mentioned above, are shown. 

It can be noted that the experimental GL parameter is almost constant at low temperatures, according to the 
mean-field description, while it shows a sharp temperature-dependence (at T jT^ « 0.8 for Tl2Ca2Ba2Cu30io, at 
T/Tc ~ 0.9 for YBa2Cu207_a, at T/Tc ~ 0.6 for MgB2) in the pre-critical region: it is in this region that the mean- 
field approximation breaks down, while our model shows its effectiveness. We would like to emphasize that the GEP 
approximation is able to reproduce the GL parameter behaviour for classes of compounds (cuprates and diborides) 
with very different microscopic structures and, probably, different 'pairing mechanism. These differences are reflected 
in the cutoff values we used to fit the data: they have the same order of magnitude (A^o = 20 for Tl2Ca2Ba2Cu30io, 
A^o = 5 for YBa2Cu207_5) for the two cuprates, while the cutoff used for the diboride (A^o = 900 for MgB2) is 
about two order of magnitude greater. In this sense, the parameter A seems to be a characteristic energy scale of the 
sample, related to its microscopic structure, and, as a consequence, to the size of the region where fluctuations cannot 
be neglected. 

Finally, we want to remark that, very close to the critical point, some universal behaviour should be expected and 
the role of thermal fluctuations becomes too important to be dealt with by the present method. In fact the GEP 
always predicts a weak first order transition at the critical point even for the neutral superfiuid (real scalar theory) |30j| . 
This is not a problem for the interpolation of the experimatal data as the difference only arises in a very narrow range 
of temperature at the transition point. Thus the method seems to be suitable for interpolating the pre-critical region 
where the order of the transition does not make any difference. Any extension to the critical range would require the 
proper inclusion of vortex fiuctuations which are expected to play a major role at criticalitv,31j. and have been recently 
shown to change the order of the transition in the GL model[23. We also mention that a variational perturbative 
extension of the gaussian approximation has been shown to be able to reach the critical regime |32| 

In summary we have shown that a polar decomposition of the complex field in the action has the merit of allowing 
the integration of the longitudinal vector field (by a further variational approximation) , is better than the Cartesian 
decomposition in and, at variance with the mean-field approximation, yields a one-parameter interpolation scheme 
which fits very well all the available data, in the whole range of the accessible temperatures. Moreover, the variational 
character of the model, allowing the study of any kind of coupling, makes the GEP a suitable tool for the inspection 
of other physical systems, such as superfiuids. 

As the GEP provides a nice way to interpolate the experimental data beyond the mean-field regime, we expect the 
method to be reliable for the description of symmetry breaking in 3 + 1 dimensions where the scalar electrodynamics 
may be regarded as a toy model for the standard electro-weak theory. 
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FIG. 1: Comparison between the polar GEP (solid line) and the Cartesian GEP dashed line. The value of the parameters are: 
m%/\% = —80, A/As = 10, cb/Xb = 0.001. Note that the origin of the <p axis coincides with the minimum of the Cartesian 
GEP to emphasize that our polar decomposition of the field leads to a lower effective potential. 




FIG. 2: The GL parameter according to Eq.lgSl for ko = 100, = 1.36 nm, Tc = 
circles are the experimental data of Ref. |26|l for Tl2Ca2Ba2Cu30io. 
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